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ABSTRACT 

Recently there has been a lot of interest in models in which gravity becomes 
strong at the TeV scale. The observed weakness of gravitational interactions 
is then explained by the existence of extra compact dimensions of space, which 
are accessible to gravity but not to Standard Model particles. In this letter 
we consider graviton emission into these extra dimensions from a hot supernova 
core. The phenomenology of SN1987A places strong constraints on this energy 
loss mechanism, allowing us to derive a bound on the fundamental Planck scale. 
For the case of two extra dimensions we obtain a very strong bound of M ^ 50 
TeV, which corresponds to a radius R ^ 0.3 /im. While there are a lot of sources 
of uncertainty associated with this bound, we find that pushing it down to the 
few- TeV range, which could in principle be probed by ground-based experiments, 
is disfavored. For three or more extra dimensions the SN1987A constraints do 
not exclude a TeV gravitational scale. 
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Recently, Arkani-Hamed, Dimopoulos and Dvali proposed a novel solution to the hierar- 
chy problem, which does not rely on either low-energy supersymmetry or compositeness of the 
Higgs boson 0, [|. They pointed out that if there exist new compact spatial dimensions, 
the fundamental (higher-dimensional) Planck scale M could be close to the electroweak scale 
Mew, thus avoiding the hierarchy. That is, gravity could become comparable in strength to 
other interactions at energies of TeV order. At distances large compared to the size of the 
compact dimensions, gravity obeys the 4-dimensional Newton's Law, with the gravitational 
constant given by 

G N = —M~ n - 2 R- n , 

where n is the number of extra dimensions, and R is their sizeQ In order to reproduce the 
measured value = 6.7 x 10 -39 GeV~ 2 with the fundamental scale M ~ 1 TeV, we require 

R ~ 2 x 10 31/n " 17 cm. (1) 

The most obvious experimental consequence of this scenario is the violation of Newton's 
Law at distances of order R. Macroscopic measurments of gravity constrain R to be less 
than about a millimeter Q. From ([!]) we see that for n = 1, R ~ 10 13 cm, so this case is 
clearly excluded. For n > 2, however, this constraint is satisfied: M = 1 TeV corresponds 
to R = 0.68 mm for n = 2 and to R = 3.0 x 10~ 12 cm for n = 6. 

Since the Standard Model provides an accurate description of the strong, electromag- 
netic and weak processes up to energies of order 100 GeV, the quarks, leptons and gauge 
bosons cannot propogate in the extra dimensions. Therefore they must be localized to a 
4-dimensional hypersurface within the full space-time. While such a localization may be 
achieved by purely field-theoretic mechanisms JIJ, the most attractive possibility is provided 
by D-branes of string theory [fj], since they naturally appear with certain degrees of freedom 
confined to them. 

If indeed gravity becomes strong at TeV scale, the higher- dimensional gravitons should 
have significant couplings to the Standard Model particles at energies accessible to current 
and near-future collider experiments |§. In particular, graviton production in e + e~ and pp 
collisions leads to events with missing energy]^, 0. The current results from searches for such 
events at LEP and Tevatron require that M be higher than 1.2 TeV for n = 2 and 610 GeV 
for n = 6 |. A number of authors have also considered the effects of the virtual graviton 
exchange on various observables |5|, |7j, as well as graviton loop effects on the anomalous 
magnetic moment of the muon M. 

Various astrophysical and cosmological processes can also be used to put constraints on 
the model. In particular, the order-of-magnitude estimates in show that the agreement of 
the observed neutrino fluxes from the supernova SN1987A with the predictions of the stellar 



collapse models |IT] requires that the fundamental scale be as high as M ~ 30 TeV for n = 2. 



Since this is one of the most stringent bounds known on M, we feel that a more detailed 



* Throughout this paper, we will assume for simplicity that the new dimensions are compactified on a 
torus of periodicity 2wR in each direction. 
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quantitative study of graviton emission from the supernova is warranted. In particular, the 
next-generation collider experiments, as well as (for the case of two extra dimensions) the 
upcoming precision measurements of gravity on short distances, could probe values of M up 
to a few TeV |4], [|, |J. Therefore it is interesting to ask whether this range is already ruled 
out by SN1987A. We will address this question in this letter. 

According to the standard theory of type-II supernovae, most of the ~ 10 53 ergs of 
gravitational binding energy released during core collapse is carried away by neutrinos. This 
hypothesis was essentially confirmed by the measurements of neutrino fluxes from SN1987A 
by Kamiokande || and 1MB |10| collaborations. These measurements allow one to put 



powerful constraints on new physics. Indeed, if there is some novel channel through which 
the core of the supernova can lose energy, the luminosity in this channel should be low 
enough to preserve the agreement of neutrino observations with theory. This idea was 



used to put the strongest experimental upper bounds on the axion mass fl2| . Here, we 
will consider the emission of the higher-dimensional gravitonsQ from the core. Once these 
particles are produced, they escape into the extra dimensions, carrying energy away with 
them. The constraint on the luminosity of this process can be converted into a bound on 
the fundamental Planck scale of the theory, M. 

During the first few seconds after collapse, the core contains neutrons, protons, electrons, 
neutrinos and thermal photons. There are a number of processes in which higher-dimensional 
gravitons can be produced. For the conditions that pertain in the core at this time (tem- 
peratures T ~ 30 — 70 MeV, densities p ~ (3 — 10) x 10 14 g cm -3 ), the dominant process is 
nucleon-nucleon "gravi-strahlung" , 

N + N ^ N + N + G, (2) 

where A" can be a neutron n or a proton p, and G is a higher-dimensional graviton. Elec- 
tromagnetic processes in which gravitons can be emitted, such as ep — > epG and ee — > eeG, 
are much less important at temperatures of interest than the strong force mediated reaction 
(@). The processes involving photons, 77 — > G, — ► eG and 7p — > pG, are also suppressed, 
since the photons have significantly smaller number density than nucleons. In this paper, 
we will only consider reaction (|[). 

To model the nucleon-nucleon interaction, we have employed the one-pion exchange ap- 
proximation. The nucleon-pion interaction Lagrangian is 



C NNn = --^d^ a N l5 ^r a N, (3) 



where N T = (p,n), r a are Pauli matrices, qa — 1.25, and F n ~ 180 MeV. While the 
validity of the one-pion exchange potential at nuclear densities is questionable, we believe it 

T It is likely that, in realistic models, the extra dimensions will contain scalar, vector, and even fermion 
fields that have substantial couplings to matter |0|. While these particles could also contribute to the 
energy loss from the core and strenghten our bound, we do not consider them here, since this contribution 
is model-dependent. 
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Figure 1: Diagrams that contribute to the gravi-strahlung process in nucleon-nucleon colli- 
sion in the one-pion exchange approximation. 



is appropriate to use it for our present purposes, since factors of order unity in the graviton 
emission rate will not have much effect on our conclusions. 

From the four-dimensional point of view, a higher- dimensional graviton appears as a tower 
of "Kaluza-Klein modes", massive particles whose mass is determined by their momentum 
in the extra dimensions and is quantized in units of 1/R. The typical energies at which 
gravitons are emitted in the core are much lower than the fundamental scale; in this limit, 
only spin-2 modes (which we will refer to as "gravitons" , h^ u ) and spin-0 modes ( "dilatons" , 
4>) couple to the Standard Model particles^]. The Lagrangian is [|15j: 



^ = -^E( rJ \ + ^). ( 4 ) 

3 

where T^ u is the conserved energy-momentum tensor of the matter, « = \Z32ttGn, oj = 
\f 3(n+2) ' anc ^ ^ ne sum * s over var i° us Kaluza-Klein (KK) modes, labelled by their mo- 
mentum in the compact dimensions j. Note that each KK mode couples to matter with 
the usual four-dimensional gravitational strength. However, for n < 4, the number of these 
modes that can be emitted at the core temperatures is very large: the mass splitting be- 
tween the modes ranges from 10~ 4 eV for n = 2 to about 10 keV for n — 4. This enormous 
multiplicity leads to enhancement of their effects, and allows us to put strong constraints on 
the model. For n > 4, the mass splitting becomes comparable to the temperatures in the 
core, so that only the first few modes can be emitted. In this case, the bound we obtain is 
very weak, and we will not discuss it here. Note also that once a particular mode is emitted, 



■fThis fact is obvious from the effective field theory point of view, developed in |Q. 
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it practically does not interact with matter. Thus, the gravitons, unlike the axions and neu- 
trinos, cannot be trapped in the core. We emphasize that the Lagrangian (f|) is completely 
independent of the details of the physics at the fundamental scale M. Therefore, our results 
will provide a model-independent bound on this scale. 

Given the Lagrangians (H) and (|§), it is a simple matter to evaluate the matrix elements 
for each of the processes in (|2]). There are 14 diagrams to evaluate, the seven "direct" 
diagrams shown in Figure [I], and seven "exchange" diagrams, with legs 3 and 4 interchanged. 
Since the temperatures in the core are of order 30-70 MeV, the nucleons are nonrelativistic. 
Naively, one would expect the diagrams (a)-(d) (and the corresponding exchange diagrams) 
to dominate in this regime. However, the leading contributions cancel due to energy and 
momentum conservation, and all the diagrams have to be included in the calculation. We 
have also made the approximation that m 2 <C 3m7vT; this is satisfied for T>6 MeV. For 
the emission of a single Kaluza-Klein graviton mode j we find 

53 \M'{nn -> nnG)\ 2 = £ \MP{pp -> ppG)\ 2 = A 1 + B 1 - 1C\ 

spin spin 

^ \M?{np -> npG)\ 2 = A? + 4B? + 4C l , (5) 

spin 

where A 3 and fv 7 are the contributions of the direct and exchange diagrams, respectively, 
and O is due to the interference between the two types of diagrams. A minus sign appears 
in the first line due to exchange of identical fermions in the final state. In np scattering, the 
contributions of exchange diagrams are enhanced because they involve charged pions, which 
couple to nucleons more strongly than neutral pions by a factor of \/2. Explicitly, we find 



\ca 4 _4 r 



1024vr m%g%G! 



0> 



45 Ft 



5127T m%g\G 



N 



45 Fi 



27 + 26 m ^ m ' + 7 m ^ m ' + l^ 2 !^ 2 ( ig + 22^4 + 4 m W X 1 

(k ■ i) 2 (k ■ iy (k ■ if \ (k ■ if (k ■ Ty , 

27 + 26 m ^ + 7 m ^ m g - 5^-5! (7 + + 2 m N m t \ 

(k-f) 2 (k-iy \k\ 2 \t\ 2 \ (k-f) 2 (k-fy) 

I { \k\ 2 \r\ 2 ^k-fjA f 22 m 2 N m 2 g m%mj" 

ik-f) 2 \k\w v {k-iy {k-iy, 



(6) 



where k = p[— p^, I = pi — Pa, and m g is the four-dimensional mass of the emitted graviton 
mode. These expressions have been averaged over the direction of the graviton. We have 
also evaluated the matrix elements for the dilaton emission, but since their contribution to 
the energy loss rate turns out to be negligible we do not present them here. 

The energy loss rate for the nucleon-nucleon gravi-strahlung is given by the following 
phase-space integral: 

e = E/ dU 1 dU 2 dU 3 dU 4 dU g S\M'\ 2 (2nyd\p 1 +p 2 ~p 3 -p 4 -p g )E g f 1 f 2 (l-f 3 )(l-f 4 ), (7) 
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where dlli = d 3 pi/ (2ii) 3 2Ei, the labels i denote the incoming (i = 1, 2) and outgoing (i = 3, 4) 
nucleons, the label g denotes the higher- dimensional graviton, S is the symmetry factor for 
identical particles in the initial and final states (5 = 1/4 for nn —>■ nn, pp —>■ pp, and S — 1 
for np — > np), and fi = [exp(Ei/T — fii/T) + are the distribution functions of the 
nucleons. Since the typical energies in our process are much higher than the splitting of the 
KK modes, we can replace the sum over these modes by an integral over the four- dimensional 
mass, according to 

£ ^ R nJ d n m = }_ QnR n J (^jM)^^ = ^M~^ G# J (m 2 )^/ 2 d m 2 , (8) 
1 

where Q n is the surface area of the unit sphere in n dimensions. 

We will assume that the nucleons in the core are non-degenerate; in this limit, fi = 
exp(yUj/T — Ei/T) and we can neglect the "blocking factors", 1 — f$ and 1 — f^. Integrating 
the matrix elements fl5|) over the phase-space in (^) and over the four- dimensional mass of 
the emitted graviton, we get 



1.7 x 10 17 erg g'V 1 [X 2 n + X 2 + 7.0X n X p ) p lA T^ eY M^ Y , n = 2 



9.4 x 10 11 erg g^V 1 [X 2 n + X 2 + 7.8X n X p ) p u T^ 5 eY M^ Y , n = 3 
5.9 x 10 6 erg g"^" 1 (x 2 n + X 2 + 8.8X n X p ) p lA T 7 M %M^, n = 4 



(9) 



where X p and X n = 1 — X p are the proton and neutron fractions in the core, px4 = 
p/(10 14 g cm" 3 ), T MeV = T/(lMeV), and M TeV = M/(lTeV). 

In order to obtain reliable bounds on the fundamental Planck scale M, one would have 
to incorporate the luminosities (|9]) into a numerical code for a protoneutron star evolution, 
and calculate the expected neutrino fluxes for various values of M. For the purposes of this 
paper, however, we will adopt a simple analytic criterion, suggested by Raffelt (see p. 504 
of [fni].) Namely, we will require that the energy loss rate to gravitons, evaluated at typical 
core conditions, does not exceed 10 19 erg g -1 s -1 . Using @ then yields, for the case of two 
extra dimensions, 

M > 0.36 TeV(X 2 + X 2 + 7.0X n X p f 25 p°f T^ 375 , n = 2. (10) 

This bound depends rather weakly on the proton fraction, core density, and the exact value of 
the maximum luminosity we allow. On the other hand, its temperature dependence is strong. 
Therefore, the main uncertainty comes from the lack of precise knowledge of temperatures in 
the core: values quoted in the literature range from 30 to 70 MeV. For a numerical estimate, 
we have assumed, conservatively, T = 30 MeV, p = 3 x 10 14 g cm -3 , and X p = 0. This 
choice of parameters yields 

M > 50 TeV, n = 2. (11) 

This corresponds to an extra dimension size of 

i? < 3 x 10" 4 mm, n = 2. (12) 
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These values of M and R are well beyond the reach of both current and near-future ground- 
based experiments. We note that for higher realistic values of temperatures and densities in 
the core the bound obtained from (|T0| ) can be as high as 250 TeV. 

In the course of our calculation, we have made a number of simplifying approximations. 
In particular, we have not considered multiple-scattering effects, which suppress emission 
of gravitons with energies below the nucleon-nucleon collision rate in the medium^ The 
deviations of the nucleon mass and nucleon-pion coupling from their vaccuum values in the 
dense medium of the core were also neglected. While at the core temperatures the nucleons 
are neither strongly degenerate nor non-degenerate, we have not included the degeneracy 
effects, such as the Pauli blocking factors for the final-state nucleons. One may also worry 
about the inadequacies of the one-pion exchange approximation at short distances, e.g. the 
effects of pion loop diagrams, other meson exchanges, and the hard core of the nucleon- 
nucleon potential. Finally, we have not considered the possibility that novel phases of nuclear 
matter (pion condensates, quark-gluon plasma, superfluidity, etc.) may form in certain 
regions of the core. 

To make a crude estimate of the uncertainties introduced by neglecting the above effects, 
we used the results for the very similar process of nucleon-nucleon axion bremstrahlung, 
which has been studied in great detail. For this process, it was shown that a simple calcu- 
lation based on the same approximations we have made here gives a good first estimate of 
the bound on axion-nucleon coupling, g a ^ 10~ 10 This is only about a factor of 3 less 



restrictive than the generally accepted bound, corresponding to about a factor of 10 uncer- 
tainty in the axion luminosity. We expect a similar result to hold for the process considered 
here, even though a much more careful study is of course necessary to confirm it. If this 
is the case, our bound on M is quite robust, since changing the graviton luminosity by a 
factor of 10 only changes the bound by 10 a25 ~ 1.8. Note that in order to push this bound 
down to the range of a few TeV, which in principle could be probed by both collider and 
gravitational experiments, our naive calculation would have to overestimate the luminosity 
by more than 3 orders of magnitude. We feel therefore that it is fair to say that this range 
is disfavored, and likely ruled out, by SN1987A. 

We have repeated our analysis for n = 3 and n = 4. With the same parameter choices 
as before, we obtain the bounds: 

M > 4 TeV, R < 4 x 10" 7 mm, n = 3; 

M > 1 TeV, R < 2 x 10~ 8 mm, n = 4. (13) 

While these bounds are quite strong, our conclusions in this case are much more optimistic. 
Given the uncertainties of our analysis, models with large extra dimensions and values of the 
fundamental scale which could be probed by near-future collider experiments are certainly 
not excluded for n > 3. 



§ This phenomenon is analogous to the Landau-Pomeranchuk-Migdal effect in the usual photon 
bremstrahlung . 
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In summary, we have found that graviton emission from SN1987A puts very strong con- 
straints on models with large extra dimensions in the case n — 2. In this case, for a 
conservative choice of the core parameters we arrive at a bound on the fundamental Planck 
scale M ^ 50 TeV, which corresponds to a radius R ^ 0.3 /im. Even though taking into 
account various uncertainties encountered in our calculation can weaken this bound, it is 
unlikely that it can be pushed down to the phenomenologically interesting range of a few 
TeV. For n = 3 and 4, we find that the fundamental scale has to be higher than about 4 and 
1 TeV, respectively. 

We are grateful to Nima Arkani-Hamed, Savas Dimopoulos, Lance Dixon, Michael Peskin, 
Georg Raffelt, Scott Thomas and Robert Wagoner for useful discussions. M. P. was supported 
by the Department of Energy under contract DE-AC03-76SF00515. S. C. was partially 
supported by NSF grant PHY-9870115. 
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